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An argument is given showing that Coulomb attraction between conduction electrons and impurity 
ions in a dilute magnetic alloy (DMA) can be disregarded, provided the system's inverse temper- 
ature f3 is replaced by an effective inverse temperature t < f3. This replacement allows to remove 
the singularity in Kondo's expression for DMA impurity resistivity and extend his theory to K. 
The extended Kondo formula agrees with experimental data on resistivity of CuFe in the range 
of low temperatures and in the neighbourhood of the resistivity minimum. 

Using an asymptotic solution of the thermodynamics of a dilute s-d system at inverse tem- 
perature t, the impurity thermodynamic functions are derived and shown to provide good agree- 
ment with experimental data on CuFe, CuCr and (LaCe)Al 2 alloys in the low-temperature range. 
The magnitude of these functions agrees with experiment and does not require rescaling as in previ- 
ous s-d theories. Nonlinear dependence of CuFe heat capacity on impurity concentration has been ac- 
counted for the first time. 
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I. INTRODUCTION 



The anomalous thermal behaviour of dilute magnetic alloys (DMA) has been the subject of extensive experimental 
and theoretical research over the past decades. The main stream of theoretical investigations (e.g. Refs. [ll^HIIllIluI 
MIM^ ITfllTlh has focused on the construction of a conductivity theory and thermodynamics of DMA on the grounds 
of the s-d Hamiltonian H s _^ introduced by Kasuya^. 

Kondo's theory of DMA impurity resistivity Ap~ has successfully explained the experimentally observed dependence 
of Ap on impurity concentration c and temperature T in the vicinity of the resistivity minimum. Unfortunately, 
the theory fails at T = 0, where Kondo's expression for Ap exhibits a logarithmic singularity. The question of removing 
this singularity is known as the "Kondo problem" . 

Substantial progress in understanding the anomalous DMA thermodynamics was made by Andrei et al££ 
who solved H s _d thermodynamics for an s-d system with indistinguishable impurities and point s-d interaction. 
Their rescaled thermodynamic functions, agree, up to a small error, with experimental impurity specific heat and mag- 
netization of (LaCc)Al 2 . The solution found in Refs. [5|j6| yields universal single-impurity curves for each thermody- 
namic function corresponding to a given value of impurity spin, which are independent of impurity concentration c. 

The shape of experimental plots of DMA thermodynamic functions, in general, varies slowly with c (e.g. Refs. [l], 
ITMTil ). meaning that their dependence on c is nonlinear. This type of dependence has been accounted for by theory 
only in exceptional cases (e.g. Refs. [HflU) . 

A different solution of H s _a thermodynamics, which treats a dilute s-d system (with a smeared s-d interaction) con- 
taining arbitrarily positioned distinguishable impurities and yields nonlinear dependence of thermodynamic functions 
on c, was presented by the author in Refs. flCilfrn . The first quantization Hamiltonian i? s _d in this approach is 

(M n \ 
H o n) ^( R « - r >) ® S *«°*i > (!) 

a=l i=l ) 

n (M) denoting the number of electrons (impurities), the antisymmetrizer with respect to electron variables 

with indices i = 1, . . . , n, and 

n J; 2 

ST* (2) 

i=l 

R Q denotes the position vector of the ath impurity, S za its spin operator and r^, a Z i stand for the respective quantities 
of the ith electron. U > represents any sufficiently regular function depending on |R Q — r<| , which allows application 

of the Feynman-Kac theorem to exp[— /3H^ This theorem was applied in Ref. [l^ to derive an upper and lower 
bound to the system's free energy per electron, 

f(H^ M \(3) := -(n/3)- 1 lnTrexp[-/3i7 s ( _" d M) ], 
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and to prove that the two bounds coaleesce in the limit of small c — > (d-lim) and, as a consequence, that 

lim d-lim f{H^ M) , /3) = lim d-lim min f(h {n ' M) (£, rj), (3) = lim f{H^ ] A {n \ 0), (3) 

n— >oc n—*oo £<n n—*oo 

where (£> ??) is thc 

mean-field Hamiltonian of an s-d system Sq with separated electron and impurity vari- 
ables. Both the electron subsystem S e and the impurity subsystem Si mp of So consist of noninteracting particles. 
According to Eq. j3]), h^ n ' M '(^,rf) is almost thermodynamically equivalent to _ff s _d in the extreme dilute limit. 
The 1-electron Hamiltonian of S c has the form 

ttf'^&v) = ^' M) + \M{e - (4) 



with 

M 

h?Mfo v) = H W g ^ _ £ E#> ® «r« (5) 

ujp denoting the multiplication operator by U{H a — rj) and r](t;) = £ — /2(C) > where 

!r(B exp(-/3/i)) 
Trexp(-/3/i) ' 



-9 /c\ / R \ Tr(i?exp(-/3fe)) 



whereas /ijm p is the 1-impurity Hamiltonian of Si mp : 

h'£ p (0=9V^S za + ~g 2 S 2 za . (7) 
The necessary condition for the minimum in Eq. © takes the for m 

€ = /3(e) (8) 

with 

M():=h(M())+H(), (9) 

M^^,, (fV do) 



pn_g(l) 



AW^Wglf"- 1 ))^, (11) 



nr™/^ 1 -* denoting the Hamiltonian of n noninteracting fermions with the 1-fermion Hamiltonian M 1 ' (cf. Ref . [l9l) . 

The mean-field thermodynamics founded on Eq. ([3]) was used in Ref. [l?] to explain the presence of thc impurity 
heat capacity peak of CuCr and (LaCe)Al 2 in the vicinity of the Kondo temperature Tk- In contradistinction to earlier 
papers (e.g. Refs. I5ll6ll20fj scaling procedures were not used and nonlinear dependence of the CuCr peak's shape on c 
was taken into account. 

One of the shortcomings of DMA theories founded on s-d type Hamiltonians is the omission of the Coulomb at- 
traction between impurity ions and conduction electrons. This problem has been treated in various ways in the past. 
Kondo introduced an additive term in his resistivity formula^ to account for these interactions. In Refs. [21I22I the equi- 
librium state of an electron gas interacting with impurity ions was studied by averaging the 1-electron Green's function 
over impurity positions. By applying this method to the 1-particle equilibrium density matrix of a quantum particle 
in a field of randomly positioned wells, representing the screened Coulomb potential at each impurity site, it was shown 
in Ref. [I?] that in the low temperature regime, a gas of such particles behaves effectively, with respect to 1-particle 
measurements, like a gas of free particles at an inverse temperature t related to the system's real temperature T 
by the equality 



t(8,T) = ^ 1 tanh((5(fc B T)" 1 ) 



(12) 
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where S = \ft,\J Muzmr 1 , U2 denoting the 2nd derivative at the well's minimum. Accordingly, the inverse tem- 
perature (3 — (fc^T) -1 of the s-d system under consideration will be subsequently replaced by t(5, T + AT) 
and S, AT will be treated as adjustable parameters. The shift AT is introduced in order to compensate omission 
of other DMA interactions in iT_d- 

In the high temperature range, t(S,T) approaches smoothly /3, whereas \imt(S, T) = 5 _1 as T — > 0. Replacement 
of /3 by t in Kondo's resistivity formula^, therefore allows to account for the Coulomb interactions between impurity 
ions and conduction electrons and to remove the singularity in his theory. In Section [TT] the resulting expression 
for DMA impurity resistivity Ap is shown to give a good fit with experimental Ap for CuFe with c = 110 ppm. 
Kondo's expression for the total DMA resistivity p is also shown to comply with experimental p under this substitution. 

The objective of subsequent sections is to study the impurity heat capacity AC, magnetization AM and susceptibil- 
ity A% of the mean-field system Sq and to adjust the constants which enter these thermodynamic functions to obtain 
the best possible agreement with experiment. For alloys with spin 1/2 impurities, such as CuFe and (LaCe)Al 2 , 
there is good agreement between theory and experiment, as regards dependence on c and T in the low-temperature 
range. In particular, nonlinear dependence of AC on iron concentration in CuFe has been accounted for successfully 
the first time. The magnitude of all thermodynamic functions agrees with experiment. For the CuCr alloys, containing 
spin 3/2 ions, agreement is slightly weaker, presumably due to the simplicity of the assumed s-d interaction in Eq. (flj). 
which permits only orbital s-wave scattering (cf. Ref. 13)- 

Computations were carried out using Wolfram's Mathematica 5.2. 



II. A SOLUTION OF THE KONDO PROBLEM 



In 1964 Kondo derived his well known formula 1 - for the impurity resistivity of DMA: 

Ap = cp A + cp M (l - 3zg 2 e F 1 InT) (13) 

where pa, PMi z are constants and z is positive for antiferromagnetic s-d interaction. Inclusion of lattice resis- 
tivity pl — aT 5 yields Kondo's expression for the total DMA resistivity/? , which provides a good fit to the c, 
T dependence of experimental data on resistivity of several DMA (cf. Ref. [0). 

The breakdown of formula (|13p at T = can be easily amended by noting that Kondo's theory takes into account 
the Coulomb attraction between conduction electrons and impurity ions simply by including an additive constant 
into the r.h.s. of Eq. (|13p . From the viewpoint of Eq. ([12")) it would be more appropriate to replace in Eq. (fT3")) 
the true inverse temperature of the alloy by t. For AT = 0, the resulting expression for Ap then takes the form 

Ap = Ap + Apim(tanh(— -)) (14) 

k B l 

and is regular at T = 0. 

Plausibility of formula (fl4|) was tested by adjusting the constants Apo, Apr, 6 to fit experimental Ap(T) data 
for CuFe, with c = 22 ppm, plotted in Ref. [23|. The function (TT4")) . for Apo = 1.455/10 9 ohm cm per ppm, Api = 
0.07/10 9 ohm cm per ppm and 8 — 2/10 4 eV, is depicted in Fig.[TJ Agreement with experiment is good. 

The experimental Ap data for CuFe with c = 22 ppm are quite typical. Resistivity measurements of a variety 
of DMA samples point to the close similarity of Ap(T)/Ap(0) curves for various alloys^ 4 -. Formula (fT4|) can be 
therefore expected to provide good agreement with experiment for a large class of DMA. 

Hamann's expression for impurity resistivity^, with t replacing (3, was also confronted with the data of Ref. |23| . 
Qualitative agreement was found. 

Performing the substitution (3 — * t in Kondo's expression for total DMA resistivity p, one obtains 

p(T) = Ap(T) + Ap 2 ftanh(^)) . (15) 

For Ap = 319.2937xl(T 4 pohm cm f Api = 2xlCT 3 pohm cm > Ap 2 = 1.065 x l(T 7 pohm cm and 8 = 2.66xl(T 4 eV, 
the function p(T) provides a good fit to the experimental plot of p(T) for CuFe with c — 1.23 x 10~ 3 (Ref. [26T ) 
and is depicted in Fig. [21 

The graphs of ksT and t(5,T)~ 1 for <5 = 2.66 x 10~ 4 eV are plotted in Fig. [3] Close similarity of the two plots 
above 15 K shows that Eq. (fT5"]) provides an extension of Kondo's formula (TT4"]) to the vicinity of K. 
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FIG. 1: Impurity resistivity Ap of CuFe, with c = 22 ppm, as given by Eq. (|14|l . for Apo = 1.455/10 9 ohm per ppm, Api = 
0.07/10 9 ohm cm per ppm and 5 = 2/10 4 eV. The points are experimental results from Ref . 12a. 
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FIG. 2: The total resistivity p of CuFe, with c = 1.23/10 3 , as given by Eq. fig]), for Ap = 319.2937/10 4 ^ohm cm, Api = 
2/10 3 pohm cm, Ap2 = 1.065/10 7 /iohm cm and 5 — 2.66/10 4 eV. The points are experimental results from Ref. |26|. 



III. MEAN-FIELD IMPURITY HEAT CAPACITY 

Expressions for the mean-field energy A[/ s -d and heat capacity AC of a dilute s-d system, relative to that of a free 
electron gas, will be derived, using the Hamiltonian h^ n,M ' (£, r/), and compared with experimental data on impu- 
rity heat capacity of CuFe 13,27 , (LaCe)Al^^ and CuCr— . The system's inverse temperature (3 will be replaced 
by t(d,T). To compensate deficiencies of -ff s _d, such as non-inclusion of the interaction between d-electrons present 
in the Anderson Hamiltonian, a shift AT > in the temperature scale will be implemented (cf. Refs. 2,28). Thus the ef- 
fective inverse temperature (EIT) 



t(S, T + AT) — 5- 1 tanh (S(k B (T + AT))- 1 ) 



(16) 
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FIG. 3: The graphs of k B T and t^T)" 1 for 5 = 2.66/10 4 eV. 



A. CuFe alloys 



The spin of the Fe ions in CuFe equals 1/2 according to Ref. 1 2 71. therefore 



/ 2 (0 = tanh (t(6, T + AT) 7 >I0 (17) 



where n = Mni, 7 = yMg, n\ denoting the number of conduction electrons per impurity. f\ defined by Eq. |TT 
is a linear function in the simplest approximation: /i(£) = bo + 61 £ + . . J^. 

The expectation energy of a spin 1/2 system Si mp containing M impurities with the 1-impurity Hamiltonian (JT]) 
equals 

tAmp = (C),,,,, = -M 2 ni£/ 2 (0 + \l 2 (18) 

imp 

and, according to Appendix B of Ref. 17, the interaction energy of electrons with the Hamiltonian he (£ ; v) equals 
AU C = (hi n ' M \Z,ri)) - (nTZllP) = -Mn/ a (0(&o + hMO + •••)> (19) 



where £ is the minimizing solution of Eq. |5]| . 
Using Eqs. Pjl. |(I5 || . ([19 ]) and the definition 



h^ M \^,r,):=^(0 + ht M H^ V ) (20) 

one obtains 

AU s . d = C/ imp + AU C + Mn£/ 2 (£) - iAfn/ 2 2 (£). (21) 

The n-electron, M-impurity spin 1/2 s-d system will be now treated as a subsystem of a sample S containing one 
mole of impurities. The energy AUs — 6.022 x 10 23 M _1 AU s .d of such sample expressed in meals, equals 

AU S = 602.2 x 38271.78A/- 1 Q7 2 - ^MV/f (£) - M 2 m/ 2 (£)(&o + 61/2 (£) +•••)). ( 22 ) 

if 60, 7, 5, £ are given in powers of eV. 
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TABLE I: 



Alloy 


c 


ni 


bo[W~ 3 VeV\ 


bi 


~/[VeV] 


M 


<5[10" 4 eV] 


AT[K] 


CuFe 


0.05% 


2000 


0.1 


-101 


0.19 


10 8 


7.14 


2.6 


CuFe 


0.1% 


1000 


0.1 


-101 


0.245 


10 9 


7.48 


3 


CuFe 


0.2% 


500 


0.1 


-101 


0.288 


10 10 


7.48 


3.9 


CuCr 


212 x 10~ 7 


10 7 /212 


1.01/m 


-631 


0.086 


36000 


lO" 13 


0.78 


CuCr 


51 x 10~ 6 


10 6 /51 


1.09/m 


-461 


0.091 


248500 


io- 13 


1.05 



o 
o 



O 

< 
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FIG. 4: Impurity heat capacity AC of CuFe given by Eq. 
The points are experimental results from Refs. IT3II27I . 



23]) . with bo, b\, 7, M, 5, AT equal to the values in Table | 



The excess heat capacity of one mole of CuFe alloy, relative to that of one mole of pure Cu, then equals 

^{dAUs dAU s d£, 



AC = 



\ dT 



d£, dT 



where £ is the unique minimizing solution of Eq. ([8]) and 



dT 



dh (dh 

' dT \ ae 



(23) 



(24) 



The mean- field AC(T + AT) curves best fitting to experimental data of Refs. Il3ll27l were obtained for the values 
of M, 7, 60, bi, S, AT given in Table U and are depicted in Fig. [4j Agreement with experiment is good, especially 
for c = 0.05%, 0.1% and below 10 K. Discrepancies at higher temperatures are presumably due to experimental error, 
which increases with temperature (e.g. Refs. I28lf29h . and to an increase of the spin values of Fe ions in this temperature 
rang o 2 ' 23 . In fact, Triplett et ali 2 ^ estimate the spin of Fe ions in CuFe to be equal 3/2. 

The minimizing solutions of Eq. ([8]) are plotted in Fig. [5l Uniqueness of these solutions is proved in the Appendix. 

It is worth emphasizing how significantly variation of c affects the shape of both experimental and theoretical 
plots of AC(T + AT) in Fig. 2J It has been suggested by some authors (e.g. Refs. [2jJ3Cj) that nonlinearity of AC 
in c observed in DMA is due to impurity-impurity interactions. The above analysis shows that this nonlinearity 
can be explained on the grounds of the s-d Hamiltonian without additional interaction terms. 

Another remarkable property of the theoretical plots of AC(T + AT) in Fig. 2] is their dependence on M. The best 
fitting values Mf of M fall in the range 1 <§C Mf <C A — 6.022 x 10 23 . The sample S can be therefore viewed as made up 
of magnetic domains, each containing Mf impurities with a favoured impurity-spin orientation, which varies, in general, 
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FIG. 5: The plots of the minimizing solutions £(T + AT) of Eq. J5J) corresponding to the AC(T + AT) curves of CuFe depicted 
in Fig. H 



from one domain to another. Experiment has confirmed existence of magnetic domains in some magnetic materi- 
als (e.g. Ref. l3ll ). The orientation of electron spins is opposite to that of impurities, as follows from Eqs. (0) , ([7j> - 

Since the solution £(T) of Eq. ((SJ decreases with decreasing T (Fig. , it follows that the ordering of impurity spins 
declines as the temperature is lowered. A similar dependence of £(T) on c can be observed. In order to find lim£(T) 
as T — > for small enough c let us recall that limt(S,T) = j3 as 5 — > 0. The graphs of £(T) for t — (3, depicted 
in Fig. 3 of Ref. [13, as well as the form of Eq. ([5]) in all cases considered therein, show that as T — ► 0, lim£(T) = 
if S = 0, AT = 0. Hence, there is no ordering of impurity spins at T = if 6 = 0, AT = 0, a picture which agrees 
with the RKKY description of interactions between impurity spins in a DMA. 



B. (LaCe)Al 2 

Bader et al<2& performed interesting measurements of AC/c on (Lai_ x Ce a: )Al2, with x — 0.0064, in external 
magnetic fields ranging from kOe to 38 kOe. For strong fields Andrei et. al., obtained good agreement of their 1- 
impurity AC (T) function with rescaled data of Ref. [28|. 

According to Refs. 14!28j, a typical Kondo effect, without any superconducting side-effects observed 
in (Lai_ a; Ce a: )Al2 with Ce content exceeding x — 0.0067. However, for x — 0.0064, Bader et. al. estimate the difference 
between the expected normal state and measured superconducting-state heat capacities as insignificant. Thus a mean- 
field normal-state theory of AC/c for (Lai-^Ce^A^, with x = 0.0064, can be expected to provide good agreement 
with experiment. 

The number of valence electrons per host atom in LaAl will be assumed 8/3. For x = 0.0064, 

0.0064 4 8 , 3742 

c = = , n-i = —c = . 

2.9936 1871' 3 3 

The ground state of the Ce ion can be described by a spin with magnitude 1 /2 and a modified Lande factor g' = 
10/7 (Ref. [3), therefore in the presence of an external magnetic field H, the system's Hamiltonian is 

H^ M \H) = - - -VbHA^^ ® *«> ( 25 ) 

a i 

I denoting the identity in L 2 (E 3 ) n and 



(26) 
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FIG. 6: The plots of AC(T + AT)/c of (Lai-^Ce^Ab in various external magnetic fields, as given by Eqs. ([23]), (l27l) 
for x — 0.0064 and bo, fei, 7, M, 8, AT, go equal to the values in Table HTl The points are experimental results from Ref. |28| . 



TABLE II: 



b [10- i VeV] 


61 


7 [10" 4 \/eV] 


M 


<5[KT 5 eV] 


AT[K] 


30 


H[kOe] 


ff 2 [10" 9 eV] 


1.9 


-13101 


3.1 


65 


3.1215 


0.39 







1.4785 


1.9 


-13101 


3.1 


64 


0.6243 


0.27 


0.006 


2 


1.50156 


1.9 


-13101 


7.285 


215 


3.1215 


0.11 


0.006 


20 


2.46843 


1.9 


-13101 


10.85 


850 


12.486 


0.1 


0.006 


38 


1.38497 



where H = goH is the effective magnetic field at each impurity site and \xb the Bohr magneton. The excess en- 
ergy AUs/c of a sample S of (LaCe)Al 2 , with x — 0.0064, expressed in joules, equals 



- l AU, 



1 , 

s = |M _1 1871 x 602.2 x 160.2x 



1 



(27) 



if bo, 7, S, £ are expressed in powers of eV. For one mole of impurities 

dAU s 8AU S d£, 



AC = 



dT 



d£ dT' 



(28) 



Adjusting the parameters bo, b\, 7, M, 8, AT, go, one obtains the best fitting to experiment AC(T + AT)/c 
curves plotted in Fig.[S] The corresponding values of the parameters are given in Table [TT1 The mean-field thermody- 
namics founded on the Hamiltonian hS n,M ' thus provides satisfactory agreement with experimental data on the field 
dependence of AC(T). 

The value of go was found by adjusting AC(T+AT)/c to experiment for H = 2 kOe, with b , b\, 7 equal to their best 
fitting values for H = and allowing only small variations of M . The smallness of go indicates the strong influence 
of the Kondo effect in the formation of a polarization cloud of conduction electrons around each impurity. The cloud 
screens each magnetic ion from interactions with other magnetic ions2£ and, as implied by the inequality go <C 1, 
also from the applied field H. 
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FIG. 7: The graph of AC(T + AT)/c of CuCr with c = 21.7 ppm as given by Eqs. 
to the values in Table [I] The points are experimental results from Ref. [29l . 



23}, JSTJ, with b , h, 7, M, <5, AT equal 



C. CuCr alloys 

Triplctt et. al. 29 performed high-precision measurements of CuCr impurity heat capacity AC for a variety of con- 
centrations. For c = 51 ppm their AC(T) peak is well defined and terminates at low temperatures with a AC jump 
which they consider to be the effect of impurity-impurity interactions. Explanation of these experimental AC(T) 
data above the jump temperature is a challenge to any s-d theory. 

According to Monod et. al^ 2 - the spin of the Cr ions in CuCr equals 3/2. Thus one finds 



and 



where 



/ 2 (0 



Ui 



7 



3e -4t 7 2 M 1 sinh (3£ 7a /^£) + sinh (t-y^H^) 



M^/ni e - 4 ^ 2M_1 cosh (3t 7 Vm£) + cosh (fry^ra^) 



imp 



_/,(M)\ 



im P / h lM) 



-M 2 n^f 2 (0 + - 7 2 



(29) 



(30) 



/4(6 = 4 7 2 



-4t 7 2 M 



1 cosh (3t7- v /ni£) 



e _ 4t7 2 M -i cosh (3t 7 ^^) + C o S h (^V^) ' 
Using Eqs. (fK)|) . (|2"Tj) . one obtains the following formula for AUs (expressed in joules) of a sample S of CuCr: 

'la..,., 1 



AC/ 5 = 602.2 x 160.2 x ( ^ + / 4 (£) 



(31) 



if b , 7, <5, £ are given in powers of eV. AC(T + AT) obtains using Eqs. 

The best fitting graphs of AC(T + AT)/c for c = 21.7 ppm and c = 51 ppm are plotted in Figs. [3 M The corre- 
sponding values of parameters 60, £>i, 7, M, 5, AT are given in Table U 

Agreement with experimental data of Ref. [29| is satisfactory, especially for c = 21.7 ppm, although not as good 
as for CuFe (Section IIII Ajl . presumably due to simplicity of the assumed s-d interaction in Eq. |T]) and variation 
of Cr spin values at higher temperatures. It has been suggested 2 -^ that for larger impurity spins the s-d interaction 
should account for the momentum dependence of s-d coupling. 
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FIG. 8: The graph of AC(T + AT)/c of CuCr with c = 51 ppm as given by Eqs. J23J, (3TJ , with bo, bi, 7, M, S, AT equal 
to the values in Table U The points are experimental results from Ref. I29I . 



TABLE III: 



6o[10- n Vc7g] 


61 


7[10" 10 ^/efg] 


M 


<5[10 _17 erg] 


AT[K] 


3o 


H[kOc] 


g 2 [10- 24 erg] 


28.8576 


-13101 


3.9237 


10 4 


2 


0.045 


0.006 


1 


15.3954 


38.4768 


-13101 


7.8474 


2 x 10 4 


3.1 


0.094 


0.006 


2 


30.7908 


48.096 


-13101 


19.6185 


2.5 x 10 4 


4 


0.18 


0.006 


5 


153.9542 


48.096 


-13101 


39.237 


3 x 10 4 


4.2 


0.23 


0.006 


10 


513.1807 



IV. MAGNETIZATION OF (LaCe)Al 2 

Various measurements of DMA impurity magnetization AM (e.g. Ref. .14 and references therein) point to a similar 
field dependence of the AM(H,T) vs. H/T curves. A typical experimental plot of AM(H,T), for (L&i- x Ce x )Al2 
with x = 0.015, can be found in Ref. [IJ. The single-impurity theoretical AM{H,T) curves found by Andrei et al££ 
fit the rescaled data of Ref. H up to a small error. 

As implied by the form of the mean-field counterpart of the Hamiltonian (|25p . viz., 

a i 

for a mole of spin 1/2 impurities 

AM = l -g'ii B (Sz a ) h w = isWtanh (t(5, T + AT){ l -g' ^ B H - j^o) (32) 

" a— 1 ~ ^ ' 

where £ is the unique solution of Eq. ijHJ) with /2(C) given by Eq. (|26|) . The resulting plots of AM(H,T) for various 
applied fields are depicted in Fig. [9l The corresponding values of 60, bi, 7, M, S, go are presented in Table HTT1 
Agreement with experiment is good in the range of low temperatures, but less satisfactory at higher T. 
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FIG. 9: The impurity magnetization [AM| of (Lai-^Ce^Ah; in various external magnetic fields as given by Eq. (I32II . with x 
0.015 and &o, bi, 7, M, 6, AT, go equal to the values in Table UlTl The points are experimental results from Ref. UA. 



MAGNETIC SUSCEPTIBILITY 



The zero-field impurity susceptibility 



Ax 



/9AM dAM <9£ 



(33) 



has been the most frequently measured property of DMA. The theory of A\, developed by Souletie et al. 15 for a DMA 
with RKKY interaction between impurities, predicts a dependence of the form 



A X (T,c) = f(T/c) 



(34) 



where / is a function independent of concentration. Felsch et al?^ have confirmed approximate validity of Eq. (|34|) 
for (Lai_ a; Ce a ;)Al2 with x ranging from 0.02 to 0.06. 

Here validity of formula (|33|) . with AM given by Eq. (|32"f . is tested on A\ experimental data for CuFe with c = 
110 ppm (Ref. 23) and (Lai-xCe^Ala with x = 0.015, 0.02 (Ref. [3). 



CuFe 



Daybell et al^ expressed their measured A\ values for CuFe in emu per gram of alloy per ppm. Formula (|32|). 
expressed in these units, takes the form 



AMcuFc = ^VsMpe tanh (t(6, T + AT)(-g' 'n B H - jV^O 



(35) 



where Mp c = 1.042032405 x 10 18 is the number of Fe ions contained in one gram of CuFe with c = 110 ppm. A possible 
fit of the resulting function A\, expressed in these units, to the experimental data on A%, for CuFe from Ref. 
is presented in Fig. [TUJ The corresponding values of the parameters are given in Table HVl Concavity of the A%(T _1 ) 
curve in Fig. [TO] appears not to be fully adjustable to concavity of the experimental plot at higher temperatures, 
however, agreement is satisfactory in the range of low temperatures. The discrepancy between theory and experiment 
at higher T is presumably due to the increase of Fe spin values with increasing T (Ref. |23| ) . 
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FIG. 10: The impurity susceptibility Ax(T _1 ) of CuFe, with c = 110 ppm, as given by Eq. (J33J) , with b , bi, 7, M, S, g' , 
AT equal to the values in Table IIVI The points are experimental results from Ref. |23 . 



TABLE IV: 



Alloy 


X 


m 


6 [10- 10 V5?g] 


61 


7 [Hr 15 yerg] 


M 


5[10 -17 erg] 


9' 


AT[K] 


CuFe 




io 5 /n 


8 


-101 


1 


10 4 


1 


1.05 


0.005 




0.015 


1592/3 


2.4048 


-13101 


392370 


200 


1 


10/7 


0.2 


(Lai_ a; Ce a: )Al2 


0.02 


1192/3 
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FIG. 11: The inverse impurity susceptibility A\{T + AT) 1 of (Lai-xCerc) Ab according to Eq. ((33} , with x, bo, bi, 7, M, S, 
g', AT equal to the values in Table HVl The points are experimental results from Ref. UA. 
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B. (LaCe)Al 2 

Felsch et alj^ performed detailed measurements of (Lai-xCe^A^ susceptibility for x ranging from 0.01 to 0.2. 
Their plots of Ax(T) -1 for x — 0.01, 0.015 are almost indistinguishable and follow a Curie- Weiss law 

Ax(T) = X cf(T)T/(T + 6), xcfT = const (36) 

for T e (0.15 K,3 K). Below 0.15 K their data on Ax(T) -1 for x = 0.015, deviate from Eq. (J3BJ) to lower values, 
contrary to measurements of A\ on AuV— and earlier theorie s 4 ' 20 ' 35 ' 36 which predict a flattening-off of Ax(T) -1 to 
higher values. For x > 0.02 the A%(T) _1 curves of Ref. [l4| no longer obey Eq. (|3"6"1) and exhibit a weak concavity. 

Eq. ®, with AM expressed by Eq. (J32]), provides a good fit to the A X (T + AT)- 1 data of Ref. [H for x = 0.01, 
0.015 and, as could be expected from Eq. ([3|), a less satisfactory adjustment for x > 0.02. The plots of A\(T + AT) -1 , 
as given by Eq. (|33|) . are depicted in Fig. [TT] for a; = 0.015, 0.02. A minor deviation from the Curie- Weiss law at very 
low temperatures, similar to the one found for AuV in Ref. l34l . can be observed. The corresponding values of bo, b\, 
7, M, S are given in Table QV] 



VI. CONCLUDING REMARKS 

The mean- field theory of dilute s-d systems presented in Ref. [H! has, in general, proved successful in providing quan- 
titative explanation of the T, c, H dependence of DMA heat capacity, magnetization and susceptibility in the range 
of low temperatures. 

It has been shown that nonlinear dependence of DMA impurity heat capacity on c can be explained in the dilute 
limit exclusively in terms of the s-d interaction, without introducing an impurity-impurity potential. 

The EIT t has improved the dependence of all thermodynamic functions on temperature and removed the singularity 
in Kondo's expression for DMA impurity resistivity. 

Deviations of the theory from experiment, in the case of heat capacity, magnetization and CuFe susceptibility, 
are presumably due to increasing spin values of impurity ions observed in some DMA at higher temperatures and 
simplicity of the s-d coupling assumed. It has been suggested that, apart from s-wave, the coupling should account for 
the d-type character of the interaction. 

Further improvement of the theory can be expected after including higher expansion terms of /i and correction terms 
to the mean-field free energy. The variation of impurity spin values at higher temperatures also suggests investigation 
of the thermodynamics of a more general s-d system containing impurity spins with various S — 1/2, 1, 3/2, . . .. 

APPENDIX 

To prove uniqueness of the minimizing solutions of the equations 

^-f(h^ M \^),P)=Q, (Al) 

|^ ( "' M) (^),/3) = 0, (A.2) 
it suffices to note that Eqs. £0}, (|X2)) take the form (cf. Ref.fiH 

C = /i«-»?) + /a(0 ) (A.3) 

V = A(C-V)- (A.4) 
Thus, rj(£) = £ — /2(C), and therefore, for r\ satisfying these equations^, 

^f(h (n ' M) (tM0),P) = |:/(^' M) (^(0),/3) =M(£-/x(/ 2 (0)-/ 2 (0)- ( A - 5 ) 
Since fi(x) = bo + b\x and by < 0, 1 /2(f) > 0, it follows that the solution £ m of Eq. ([8]) is unique and 

^f(h {nM HH,r 1 (O),f3)<0 for£<£ m , 
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